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Abstract 

The classical drift diffusion (DD) model of spin transport treats spin relaxation via an empirical 
parameter known as the "spin diffusion length". According to this model, the ensemble averaged 
spin of electrons drifting and diffusing in a solid decays exponentially with distance due to spin 
dephasing interactions. The characteristic length scale associated with this decay is the spin 
diffusion length. The DD model also predicts that this length is different for "upstream" electrons 
traveling in a decelerating electric field than for "downstream" electrons traveling in an accelerating 
field. However this picture ignores energy quantization in confined systems (e.g. quantum wires) 
and therefore fails to capture the non-trivial influence of subband structure on spin relaxation. Here 
we highlight this influence by simulating upstream spin transport in a multi-subband quantum wire, 
in the presence of D'yakonov-PereF spin relaxation, using a semi-classical model that accounts for 
the subband structure rigorously. We find that upstream spin transport has a complex dynamics 
that defies the simplistic definition of a "spin diffusion length". In fact, spin does not decay 
exponentially or even monotonically with distance, and the drift diffusion picture fails to explain 
the qualitative behavior, let alone predict quantitative features accurately. Unrelated to spin 
transport, we also find that upstream electrons undergo a "population inversion" as a consequence 
of the energy dependence of the density of states in a quasi one-dimensional structure. 

PACS numbers: 72.25.Dc, 85.75.Hh, 73.21.Hb, 85.35.Ds 
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I. INTRODUCTION 



Spin transport in semiconductor structures is a subject of much interest from the perspec- 
tive of both fundamental physics and device applications. A number of different formalisms 



have been used to study this problem, primary among which are a classical drift diffusion 
approach ll Ij, |3J] , a kinetic theory approach |4j, and a microscopic semiclassical approach 
J, Q]. The central result of the drift diffusion approach is a differential 
equation that describes the spatial and temporal evolution of carriers with a certain spin 
polarization n a . Ref. JjJ derived this equation for a number of special cases starting from 
the Wigner distribution function. In a coordinate system where the x-axis coincides with 
the direction of electric field driving transport, this equation is of the form: 
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D is the diffusion coefficient, and A and B are dyadics (9-component tensors) that depend 
on D, the mobility \i and the spin orbit interaction strength in the material. 

Solutions of Equation with appropriate boundary conditions, predict that the ensem- 



ble averaged spin |(S) 
with x according to: 

where 



a/ (S x ) 2 (x) + (S y ) 2 (x) + (S z ) 2 (x) should decay exponentially 
|(S)|(x) = |(S}|(0)e-^ (4) 



1 
I 



2D 



+ 



(fi) + C2 



(5) 



Here E is the strength of the driving electric field and C is a parameter related to the spin 
orbit interaction strength. 

The quantity L is the characteristic length over which | (S) | decays to 1/e times its original 
value. Therefore, it is defined as the "spin diffusion length". Equation (jSJ) clearly shows that 



spin diffusion length depends on the sign of the electric field E. It is smaller for upstream 
transport (when E is positive) than for downstream transport (when E is negative). 

This difference assumes importance in the context of spin injection from a metallic fer- 
romagnet into a semiconducting paramagnet. Ref. ]J pointed out that the spin injection 
efficiency across the interface between these materials depends on the difference between the 
quantities L s /a s and L m /a m where L s is the spin diffusion length in the semiconductor, a s 
is the conductivity of the semiconductor, a m is the conductivity of the metallic ferromagnet, 
and L m is the spin diffusion length in the metallic ferromagnet. Generally, o m >> a s . How- 
ever, at sufficiently high retarding electric field, L s « L m , so that L s /a s ~ L m /a m . When 
this e q nahty is estabhshed, the S p m injection efficiency is m axnn 1Z ed. Thus, ref. B, claimed 
that it is possible to circumvent the infamous "conductivity mismatch" problem j 1 3lj . which 
inhibits efficient spin injection across a metal- semiconductor interface, by applying a high 
retarding electric field in the semiconductor close to the interface. A tunnel barrier between 
the f_et and sendeondnctor Q „ r a Schott ky hazier QQ at the interface does 
essentially this and therefore improves spin injection. 

The result of ref. J]J depends on the validity of the drift diffusion model and Equation (jlj 
which predicts an exponential decay of spin polarization in space. Without the exponential 
decay, one cannot even define a "spin diffusion length" L. The question then is whether 
one expects to see the exponential decay under all circumstances, particularly in quantum 
confined structures such as quantum wires. The answer to this question is in the negative. 
Equation {I}, and similar equations derived within the drift diffusion model, do not account 
for energy quantization in quantum confined systems and neglect the influence of subband 
structure on spin depolarization. This is a serious shortcoming since in a semiconductor 
quantum wire, the spin orbit interaction strength is different in different subbands. It is 
this difference that results in D'yakonov-Perel' (D-P) spin relaxation in quantum wires. 
Without this difference, the D-P relaxation will be completely absent in quantum wires and 
the corresponding spin diffusion length will be always infinite The suband structure is 
therefore vital to spin relaxation. 



II. SEMICLASSICAL MODEL OF SPIN RELAXATION 



In this paper, we have studied spin relaxation using a microscopic semiclassical model 
that is derived from the Liouville equation for the spin density matrix Q, 0|. Our model 
has been described in ref. |10| and wil not be repeated here. This model allows us to study 
D'yakonov-Perel' spin relaxation taking into account the detailed subband structure in the 
system being studied. 

In technologically important semiconductors, such as GaAs, spin relaxation is dominated 
by the D'yakonov-Perel' (D-P) mechanism 12]. This mechanism arises from the Dresselhaus 
19j and Rashba |20J spin orbit interactions that act as momentum dependent effective mag- 
netic fields B(k). An electron's spin polarization vector S precesses about B(k) according 
to the equation 

f = n(k) x s (6) 

where Q(k) is the angular frequency of spin precession and is related to B(k) as f2(k) = 
(e/m*)B(k), where m* is the electron's effective mass. If the direction of B(k) changes 
randomly due to carrier scattering which changes k, then ensemble averaging over the spins 
of a large number of electrons will lead to a decay of the ensemble averaged spin in space and 
time. This is the physics of the D-P relaxation in bulk and quantum wells. In a quantum 
wire, the direction of k never changes (it is always along the axis of the wire) in spite of 
scattering. Nevertheless, there can be D-P relaxation in a multi- subband quantum wire, as 
we explain in the next paragraphs. 

We will consider a quantum wire of rectangular cross-section with its axis along the [100] 
crystallographic orientation (which we label the x-axis), and a symmetry breaking electric 
field E y is applied along the y-axis to induce the Rashba interaction (refer Fig. [TJ. Then, 
the components of the vector f2(k) due to the Dresselhaus and Rashba interactions are given 
by 
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where 042 and a 4 6 are material constants, (m, n) are the transverse subband indices, k x is the 
wavevector along the axis of the quantum wire, and W z , W y are the transverse dimensions 
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Thus, B (k) lies in the x-z plane and subtends an angle 9 with the wire axis (x-axis) 
given by 



Note from the above that in any given subband in a quantum wire, the direction of B(k) 
is fixed, irrespective of the magnitude of the wavevector k x , since 9 is independent of k x . As 
a result, there is no D-P relaxation in any given subband, even in the presence of scattering. 

However, 9 is different in different subbands because the Dresselhaus interaction is differ- 
ent in different subbands. Consequently, as electrons transition between subbands because 
of inter-subband scattering, the angle 9, and therefore the direction of the effective magnetic 
field B(k), changes. This causes D-P relaxation in a multi- subband quantum wire. Since 
spins precess about different axes in different subbands, ensemble averaging over electrons 
in all subbands results in a gradual decay of the net spin polarization. Thus, there is no D-P 
spin relaxation in a quantum wire if a single subband is occupied, but it is present if multiple 
subbands are occupied and inter-subband scattering occurs. This was shown rigorously in 
ref. Q- 

The subband structure is therefore critical to D-P spin relaxation in a quantum wire. In 
fact, if a situation arises whereby all electrons transition to a single subband and remain 
there, further spin relaxation due to the D-P mechanism will cease thereafter. In this case, 
spin no longer decays, let alone decay exponentially with distance. Hence, spin depolariza- 
tion (or spin relaxation) cannot be parameterized by a constant spin diffusion length. 

The rest of this paper is organized as follows. In the next section, we describe our model 
system, followed by results and discussions in section IV. Finally, we conclude in section V. 

III. MODEL OF UPSTREAM SPIN TRANSPORT 

We consider a non centro-symmetric (e.g. GaAs) quantum wire with axis along [100] 
crystallographic direction. We choose a three dimensional Cartesian coordinate system with 



9 = arctan 
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x coinciding with the axis of the quantum wire (refer Fig. The structure is of length 
L x = 1.005 fim with rectangular cross section: W y = 4 nm and W z = 30 nm. A metal gate 
is placed on the top (not shown in Fig. P) to induce the symmetry breaking electric field 
E y y, which causes the Rashba interaction. In a quantum wire defined by split Schottky 
gates on a two-dimensional electron gas, E y y arises naturally because of the triangular 
potential confining carriers near the heterointerface. We assume E y = 100 kV/cm (lOj). In 
addition, there is another electric field —E x x (E x > 0) which drives transport along the 
axis of the quantum wire. Consider the case when spin polarized monochromatic electrons 
are constantly injected into the channel at x = Xq = 1 \xm with injection velocities along 
—x. If these electrons occupy only the lowest subband at all times, then there will be no 
D'yakonov-Perel' relaxation [171 ] . Therefore, in order to study multisubband effect on spin 
dephasing of upstream electrons, we inject them with enough energy (Eq) that they initially 
occupy multiple subbands. We ignore any thermal broadening of injection energy [3] since 
Eq » ksT for the range of temperature (T) considered, k& being Boltzmann constant. Let 
Ei denote the energy at the bottom of ith subband (i = 1, 2, • •• n, n+1, ■■ ■ etc.). We place E 
between the n-th and (n + l)-th subband bottoms as shown in Fig. |21 In other words, E n < 
E < E n+ i. We assume that the injected electrons, each with energy E , are distributed 
uniformly over the lowest n subbands. In other wordes, at time t = 0, electron population 
of the i th (1 < i < n) subband is given by Ni(x,t = 0) = (N /n) S(x — x ) 6(E — E ) 
where iVo is the total number of injected electrons and E denotes their energies. At any 
subsequent time t, these distributions spread out in space (x < x ), as well as in energy, due 
to interaction of the injected electrons with the electric field E x and numerous scattering 
events. Relative population of electrons among different subbands will change as well due to 
intersubband scattering events. Upstream electrons originally injected into, say, subband i 
with velocity —ViX (vi > 0), gradually slow down because of scattering and the decelerating 
electric field. They change their direction of motion (i.e. become downstream) beyond a 
distance \x~l — Xq\ measured from the injection point Xq. Thus, no electrons will be found 
in the i-th subband beyond x~i. Note that the value of \~x~l — xq\ depends on three factors: 
the initial injection velocity into subband i, the decelerating electric field and the scattering 
history. On the other hand, the 'classical turning distance' of monochromatic electrons 



injected into the i-ih subband with energy E for a given electric field E x is given by 



where Ei is the energy at the bottom of the i-th subband and V{ is the injection velocity in 
the z-th subband. Note that does not depend on scattering history and Xi = xi in ballistic 
transport. 

Clearly v n < t> n _i < ■ ■ ■ < v 2 < v i for a given E (see Fig. |2J. Thus, for a given 
channel electric field E x , x n — xq =min{xj — xo},i = 1,2, ...,n. Hence we concentrate 
on the region (x n ,xo) where, almost all injected electrons are upstream electrons. In the 
simulation, velocity of every electron is tracked and as soon as an electron alters direction 
and goes downstream (i.e. its velocity becomes positive) it is ignored by the simulator 
and another upstream electron is simultaneously injected from x = Xq randomly in any of 
the n lowest subbands with equal probability. This process is continued for a sufficiently 
long time till electron distributions over different subbands, Ni(x,t),i = 1,2, ...n, no longer 
change with time. Under this condition we say that steady state is achieved for the upstream 
electrons. This steady state electron distribution is extended from x = x n to x = xq and 
heavily skewed near the region x = Xq. This steady state distribution of upstream electrons 
does not represent the local equilibrium electron distribution because of two reasons: (a) 
upstream electrons are constantly injected into the channel at x — Xq, this is the reason why 
the distribution is skewed near x = xq and (b) we exclude any downstream electron from the 
distribution. At local equilibrium, there will be of course both upstream and downstream 
electrons in the distribution. 

The model above allows us to separate upstream electrons from downstream electrons 
and therefore permits us to study upstream electrons in isolation. Of course, in a real 
quantum wire, both upstream and downstream electrons will be present at any time, even 
in the presence of a strong electric field, since there will be always some non-vanishing 
contribution of back-scattered electrons to the upstream population. 

The semiclassical model and the simulator used to simulate spin transport have been 
described in ref. |l0|. Based on that model, at steady state, the magnitude of the ensemble 
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averaged spin vector at any position x inside the channel is given by 



E^ix) (S ix ) (x)\ + J2Ni(x) (S iy ) (x) + J2Ni(x) (S iz ) (x) 
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Here (5*^) (x), ( = x, y, z, denotes the ensemble average of ( component of spin at position 
x. Subscript i implies that ensemble averaging is carried out over electrons only in the % th 
subband. The above equation can be simplified to 

11 \ 2 n 

\(S t )(x)\ - 2 £ NiWNjix) |(S,)(x)| |(S,)(x)| sin 2 ^ 
^ (12) 

where (Sj) (x) = (Si x ) (x)x + (Si y ) (x)y + (Si z ) (x)z and 9ij(x) is the angle between (Si) (x) 
and (Sj) (x). Note that in absence of any intersubband scattering event, |(S;)| (x) = 1 for 
all x (i.e. initial spin polarization of the injected electrons) Q- Simulation results that we 
present in the next section can be understood using Equation (IT2j) . 



IV. RESULTS AND DISCUSSION 



We examine how ensemble averaged spin polarization of upstream electrons | (S) | (x) 
varies in space for different values of driving electric field E x and injection energy (E ) for 
a fixed lattice temperature T. We vary E x in the range 0.5 — 2kV/cm for constant injection 
energy 426 meV and T = 30 K, where Eq is measured from the bulk conduction band energy 
as shown in Fig. El The lowest subband bottom is 351 meV above the bulk conduction 
band edge. We also present results corresponding to Eq = 441 meV with E x = lkV/cm 
and T = 30K. In all cases mentioned above, injection energies lie between subband 3 and 
subband 4. Injected electrons are equally distributed among the three lowest subbands 
initially. Obviously, this corresponds to a non-equilibrium situation. All injected electrons 
are 100% spin polarized transverse to the wire axis (i.e either y or z). 

Figures OKI EJ and [TBI show how ensemble averaged spin components (S x )(x), (S y )(x), 
(S z )(x) and |(S)|(x) of upstream electrons evolve over space. Figures OH show the influence 
of the driving electric field on spin relaxation, Fig. ITT1 shows the influence of initial injection 



energy and Fig. ITBI shows the influence of the intial spin polarization. It is evident that nei- 
ther the driving electric field, nor the initial injection energy, nor the intial spin polarization 
has any significant effect on spin relaxation. Note that |(S)|(x) does not decay exponen- 
tially with distance, contrary to Equation Spatial distribution of electrons over different 
subbands is shown in Fig. UJ - and Fig. The classical turning point of electrons in 
the third subband (£3) has been indicated in each case. Fig. 171 - ITU1 show the influence of 
the driving electric field and Fig. shows the influence of initial injection energy on the 
spatial evolution of subband population. As expected, \xo — x%\ decreases with increasing 
electric field in accordance with Equation (|10p. Note that at low electric field (Figs. 7 and 
8) X3 ~ £3 since all subbands are getting nearly depopulated of "upstream" electrons at 
x = £3. Recall that X3 = X3 only if transport is ballistic; therefore we can conclude that 
upstream transport is nearly ballistic in the range |xo — x^\ when E x < 1 kV/cm. At high 
electric field, when E x > 1.5 kV/cm, (Fig. 10) {x^ — x \ > \x 3 — x \. This indicates that there 
are many upstream electrons even beyond the classical turning point. It can only happen 
if there is significant scattering that drives electrons against the electric field, making them 
go beyond the classical turning point. We can also deduce that most of these scattering 
events imparts momentum to the carriers to aid upstream motion rather than oppose it, 
since |xi — xo| > [x^ — xq\- This behavior is a consequence of the precise nature of the 
scattering events and would not have been accessible in drift-diffusion models that typically 
treat scattering via a relaxation time approximation. 

A. Population inversion of upstream electrons 

Note that even though electrons are injected equally into all three subbands, most elec- 
trons end up in subband 3 - the highest subband occupied initially - soon after injection. 
Beyond a certain distance (x = x sca t ~ 0.9 fim) subbands 1 and 2 become virtually de- 
populated. This feature is very counter-intuitive and represents a population inversion of 
upstream electrons! It can be understood as follows: scattering rate of an electron with en- 
ergy E is proportional to the density of the final state. In a quantum wire, density of states 
has 1 / \[E — ~E~i dependence where Ei is the energy at the bottom of the ith subband. As the 
injected electrons move upstream they gradually cool down and their energies approach the 
energy at the bottom of subband 3 (E3). To visualize this, imagine the horizontal line Eq in 



Fig. |21 sliding down with passage of time. As E 3 is approached, electrons will increasingly 
scatter into subband 3 since the density of final state in subband 3 is increasing rapidly. To 
scatter into a final state in subband 2 or 1 that has the same density of state as in subband 
3 will require a much larger change in energy and hence a much more energetic phonon 
which is rare since the phonons obey Bose Einstein statistics. Therefore, subband 3 is the 
overwhelmingly preferred destination and this preference increases rapidly as electrons cool 
further. Consequently, beyond a certain distance, virtually all electrons are scattered to 
subband 3 leaving subbands 1 and 2 depleted. This feature is a peculiarity of quasi one- 
dimensional system and will not be observed in bulk or quantum wells. Exact values of x% 
and x scat depend on injection energy and electric field. In the field range 0.5 — 1.5kV/cm 
and injection energy 426 meV, \xs\ > \x sca t\. However, for higher values of electric field (e.g. 
2kV/cm) or smaller values of injection energies, electrons reach classical turning point even 
before subbands 1 and 2 get depopulated. 

Because of electron bunching in subband 3, spin dephasing in the region (x sca t, x ) is 
governed by Equation (6) with n — 3. We observe a few subdued oscillations in |(S)| (x) in 
this region because of the "sine term" in Equation (|12p. However, in the region (xs,x scat ) 
subbands 1 and 2 are almost depopulated. Therefore, there is no D-P relaxation in the 
interval (x3,x scat ) since only a single subband is occupied jl^. Consequently, the ensemble 
averaged spin assumes a constant value |(S 3 )| < 1 and does not change any more. Thus 
in this region, one can say that spin dephasing length becomes infinite. It should be noted 
that it is meaningless to study spin dephasing in the region x < x 3 because electrons do not 
even reach this region. 



V. CONCLUSION 



In this paper, we have used a semiclassical model to study spin dephasing of upstream 
electrons in a quantum wire, taking into account the subband formation. We showed that 
the subband structure gives rise to rich features in the spin dephasing characteristics of 
upstream electrons that cannot be captured in models which fail to account for the precise 
physics of spin dephasing and the fact that it is different in different subbands. Because 
spin relaxation in a multi-subband quantum wire is non-exponential (even non-monotonic) 
in space, it does not make sense to invoke a "spin diffusion length", let alone use such a 



heuristic parameter to model spin dephasing. 

Finally, we have found a population inversion effect for upstream electrons. It is possible 
that downstream electrons also experience a similar population inversion. This scenario is 
currently being investigated. 
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Figure captions 



Figure 1. A quantum wire structure of length L = 1.005/xm with rectangular cross 
section 30 nm x 4 nm. A top gate (not shown) applies a symmetry breaking electric 
field E y to induce the Rashba interaction. A battery (not shown) applies an electric field 
—E x x (E x > 0), along the channel. Monochromatic spin polarized electrons are injected 
at x = xo = lfim with injection velocity — tw/. These electrons travel along — x {upstream 
electrons) until their direction of motion is reversed due to the electric field —E x x. We 
investigate spin dephasing of these upstream electrons. 

Figure 2. Subband energy dispersion in the quantum wire. 

Figure 3. Spatial variation of ensemble averaged spin components for driving electric field 
E x = 0.5kV/cm at steady state. Lattice temperature is 30 K, injection energy E = 426 
meV. Electrons are injected with equal probability into the three lowest subbands. Classical 
turning point of subband 3 electrons is denoted by x% and x sca t indicates the point along 
the channel axis where subbands 1 and 2 gets virtually depopulated. Injected electrons are 
y polarized and x = x = 1 jim is the point of injection. 

Figure 4. Spatial variation of ensemble averaged spin components for driving electric 
field E x = lkV/cm at steady state. Other conditions are same as in Figure 3. 

Figure 5. Spatial variation of ensemble averaged spin components for driving electric 
field E x = 1.5kV/cm at steady state. Other conditions are same as in Figure 3. 

Figure 6. Spatial variation of ensemble averaged spin components for driving electric 
field E x = 2kV/cm at steady state. Other conditions are same as in Figure 3. 

Figure 7. Spatial variation of electron population over different subbands at steady state 
for driving electric field E x = 0.5 kV/cm. Other conditions are same as before. 



Figure 8. Spatial variation of electron population over different subbands at steady state 
for driving electric field E x — 1 kV/cm. Other conditions are same as before. 

Figure 9. Spatial variation of electron population over different subbands at steady state 
for driving electric field E x = 1.5 kV/cm. Other conditions are same as before. 

Figure 10. Spatial variation of electron population in different subbands at steady state 
for driving electric field E x — 2 kV/cm. Other conditions are same as before. 

Figure 11. Spatial variation of ensemble averaged spin components for E = 441meV, 
E x = lkV/cm and lattice temperature T = 30K. Injected electrons are y polarized. 

Figure 12. Spatial variation of electron population over different subbands at steady 
state for E = 441meV, E x = lkV/cm and lattice temperature T = 30K. Injected electrons 
are y polarized. 

Figure 13. Spatial variation of ensemble averaged spin components for Eq = 426meV, 
E x = lkV/cm and lattice temperature T = 30K. Injected electrons are z polarized. 
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